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Abstract 

A brief  review  of  the  method  of  measuring  residual  stresses  in  poly- 


crystalline materials  with  X-rays  is  given.  The  effect  of  counting  statistics 
on  precision  is  discussed  as  well  as  factors  that  affect  accuracy.  Beam 
penetration,  stress  gradients,  and  the  form  of  the  stress  tensor  can  each 
seriously  affect  accuracy,  if  traditional  methods  are  employed,  and  practical 
procedures  to  deal  with  such  situations  are  outlined. 


°J.  B.  Cohen  is  Frank  C.  Engelhart  Professor  of  Materials  Science  and 
Engineering  and  H.  Dolle  is  a Post-Doctoral  Visiting  Scholar. 
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1.  Introduction 

Knowledge  of  stresses  (applied  or  residual)  is  quite  ing>ortant  in  an 
industrial  environment.  Crack  initiation  and  propagation  in  static  or  fatigue  load- 
ing, or  in  stress  corrosion,  can  be  impeded  by  compressive  stresses  normal  to 
the  crack,  and  greatly  accelerated  by  tensile  stresses.  Residual  stresses  * 

can  be  produced  during  quenching  (because  of  the  different  cooling  rates 
of  the  surface  and  the  Interior  and  any  volume  changes  due  to  transformations), 
during  manufacturing  (due  to  the  difference  in  deformation  of  the  surface  and 
the  interior  - as  in  straightening,  machining,  rolling  or  shot  peenlng),  or 
due  to  the  different  response  of  particles  and  matrix  in  a multiphase  material. 

During  welding,  the  resistance  of  the  cold  base  metal  to  the  solidifying  (and 
shrinking)  weld  pool  will  produce  stresses  in  both  regions.  These  stress  pat- 
terns can  and  do  change  in  service,  and  there  is  increased  interest  in  monitoring 
these.  For  a recent  overview,  see  ref.  1. 

Apart  from  deliberately  induced  residual  stresses,  the  measured  stresses 
can  develop  due  to  applied  loads,  or  even  in  a quasi-uniform  stress  field,  from 
differences  in  response  of  the  surface  and  the  interior,  one  grain  with  respect 
to  another,  or  one  phase  with  respect  to  another.  In  the  case  of  residual 
stresses  the  tensile  and  compressive  cong>onents  must  balance  across  a section, 
but  to  see  this  it  may  be  necessary  to  examine  values  in  different  phases. 

This  is  possible  with  X-rays,  but  It  has  not  often  been  of  sufficient  Interest; 
only  the  stresses  In  the  major  constituent  are  usually  examined.  Also,  even  in 
a single  phase  material  there  may  be  situations  where  it  is  not  possible  to 
remove  material  to  examine  gradients  because  of  the  geometry  of  the  piece. 

Actually,  a variety  of  non-destructive  techniques,  based  on  the  acoustic 
and  magnetic  response  of  a material  to  stresses,  are  being  evaluated  for  this 
purpose,  but  so  far  they  have  proved  to  be  too  sensitive  to  microstructural 


(2) 

variations  such  as  changes  in  texture  . The  X-ray  method,  which  was  first 

. (3) 

proposed  in  the  1920  s by  Lester  and  Aborn  , has  withstood  the  test  of  time 
in  industry;  it  is  the  subject  of  this  review,  in  ^lich  we  shall  examine  not 
only  the  method,  but  its  limitations,  accuracy  and  precision,  as  well  as  the 
latest  developments  which  have  made  it  of  considerable  use  in  field  situations. 
While  there  are  a number  of  ways  of  examining  strains  with  X-rays,  we  shall 
concentrate  only  on  those  that  are  widely  employed  in  industry.  To  prove  their 
validity,  these  approaches  have  been  compared  to  measurements  of  distortion 
after  mechanical  dissection^. 

The  basic  principle  that  is  employed  is  quite  simple:  An  interplanar 
spacing  (d)  in  the  material  serves  as  an  internal  strain  gauge.  As  shown  in 
Fig.  la,  in  a polycrystalline  specimen  in  the  normal  position  on  a diffractom- 
eter, only  those  grains  with  planes  nearly  parallel  to  the  surface  form  a dif- 
fraction peak.  If  there  are  compressive  stresses  in  the  near-surface  regions, 
the  "d"  spacing  is  increased  for  these  planes  due  to  the  effect  of  Poisson's 
ratio.  When  the  specimen  is  tilted  with  respect  to  the  X-ray  beam  (or  the  beam 
to  the  specimen),  Fig.  lb,  the  diffracting  grains  have  their  planes  more  nearly 
perpendicular  to  the  stresses  in  the  surface  and  their  spacing  is  less,  or 
even  decreased.  (Components  of  the  surface  stress  parallel  and  perpendicular  to 
the  planes  now  affect  this  spacing.)  The  26  positions  of  diffracting  peaks  from 
such  planes  and  Bragg's  law  are  all  that  are  required  to  examine  the  strain 
tensor,  from  which  the  surface  stresses  may  be  determined. 

2.  The  Basic  Equation 

In  what  follows,  we  shall  employ  the  axial  system,  L^,  for  the  measurements 

(a  "laboratory"  system)  and  axes  JP^for  the  specimen,  as  shown  in  Fig.  2.  Primed 

quantities  will  refer  to  measurements  in  the  laboratory  system,  unprimed  in 

the  sample  system.  With  d the  interplanar  spacing  at  tilts  and  d the 

o 
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value  for  the  strain  free  material,  the  strain  along  can  be  written  : 

^ 2 2 2 2 2 
e33  " e^jCos  «sin  \|r  + e12sin2$sin  ^ + e13cos$sin2\(t  + e22sin  $stn  t 

o 

2 


+ e23sin$sin2^  + e-^cos  ^ 


33 


(1) 


Once  the  strains  are  determined  (by  the  procedures  to  be  described  below). 


the  stresses,  a^j»  can  be  calculated  from: 


[e, j “ 6 


s^kx) 


(e-i-t  + e + s ii) ] • 


ij  %S2(hk£)  L6ij  °ij  %S2(hk£)S1(hki)  Vfcll  T 622  33 


(2) 


The  term  6^  is  the  Kronecker  delta  function,  and  since  ■ 0,  i + j; 
the  shear  strains  alone  determine  the  shear  stresses. 

The  St  are  X-ray  elastic  constants  for  a particular  hk£  reflection. 

For  an  elastically  isotropic  solid,  1;  S2(hk£)  - , S^(hk/)  - (-^),  where 

y,  is  Poisson's  ratio  and  E is  Young's  modulus.  Theory  exists  to  calculate  the 
constants  (for  each  hk£  reflection)  for  an  anisotropic  crystal  coupled  to  an 
isotropic  matrix,  and  this  theory  shows  that  the  values  are  within  a few  per- 
cent of  the  average  of  the  well  known  Reuss  and  Voight  formulae,  which  assume, 
respectively,  constant  average  stress  and  constant  average  strain  for  all 
grains  in  the  material.  In  addition,  theory  is  developing  for  an  elastically 
anisotropic  material  with  severe  texture.  (This  entire  area  has  recently  been 
reviewed  by  one  of  us.  Ref.  5.)  In  the  latter  case,  the  elastic  constants  associ- 
ated with  a particular  hki  reflection  may  vary  with  $ and  | tilt.  It  has  also 
been  well  established  that  the  X-ray  elastic  constants  can  vary  appreciably 
with  amount  of  plastic  deformation;  changes  as  large  as  20-40  pet.  have  been 
reported.  This  phenomenon  is  not  understood  at  the  moment  and  can  obviously 
severely  affect  the  accuracy  of  any  stress  measurement  if  calculated  values 
are  assumed  to  be  independent  of  prior  strain.  Fortunately,  it  is  possible  to 
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test  for  this  problem,  and,  in  many  cases,  to  measure  the  effective  constants. 
We  discuss  this  further  below.  In  addition,  it  is  often  the  relative  change  in 
stress  with  some  change  in  condition  that  is  important,  in  which  case  the  con- 
stants employed  are  not  as  Important. 

We  turn  now  to  the  various  methods,  after  which  we  will  discuss  other 
errors  that  affect  accuracy,  and  also  precision. 

■'  3.  Surface  Stresses  Only 

As  a free  surface  cannot  support  a stress  normal  to  it  (o^  " “ o23  " 

0),  the  stress  tensor  in  the  sample  system  shown  in  Fig.  2 is: 

/ a -.i  a,~  0 


yO  0 0 ] (3a) 

Of  course,  this  tensor  may  be  written  in  terms  of  principal  axes  in  a rotated 
axial  system: 

/ 

0 


(3b) 


Then,  it  can  be  shown  that: 

®33  “ -^5 — ~ “ % S2(hkjO[0]cos2$  + o2sin2$]sin2\|f  + S1(hkX)[®1  + c2].  (4) 

o 

o 

Note  that  at.  * ■ 0 , a comparison  of  d.  to  d immediately  indicates  whether 

T $\| f O 

the  sum  of  the  principal  stresses  is  positive  or  negative. 

Subtracting  the  value  for  ft  * 0°: 
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The  stress  in  any  direction,  $,  on  the  surface  can  be  determined  simply 

2 

from  the  slope  of  d vs  sin  The  stress-free  value,  dQ,  is  not  needed.  A 
typical  plot  is  shown  in  Fig.  3.  The  principal  stresses  and  their  directions 
can  be  obtained although  their  $ directions  are  often  known  from  the  nature 
of  the  applied  stress  system.  By  application  of  known  elastic  stresses  on  an 
identical  specimen,  the  effective  elastic  constants  can  be  obtained  from  Eqns.  4 
or  5. 

2 

It  cannot  be  emphasized  too  strongly  that  the  linearity  of  "d"  vs.  sin  ^ 
must  be  confirmed  in  each  situation,  but  if  this  is  the  case,  the  equation  can 
be  further  simplified,  with  the  aid  of  Bragg's  law,  to  involve  only  two  tilts, 
f • 0°  and,  typically,  ■ 45°  or  60°: 


a.  “ KA20  (in  degrees), 
9 


with: 


K“  180  S2(hkf)  sin2^  cot^(%-0  + (6b) 

In  this  "two  exposure  method,"  a high  angle  peak  is  chosen,  and  the  ^ tilt 

is  made  as  large  as  possible,  to  minimize  K and  hence  maximize  the  sensitivity 

of  the  peak  shift  to  the  stress.  For  example,  for  steel,  CrK  radiation,  the 

al 

211  reflection  at  « 156°20,  and  a ^-tilt  of  45°,  a + O.l°20  shift  corresponds 

to  a stress  of  « -60  MPa  (-8600  psi).  Shifts  as  small  as  0.01  - O.O2°20  can  be 

2 

detected.  If  it  has  been  verified  that  d vs  sin  ^ is  indeed  linear,  it  is 
possible  to  develop  simple  portable  equipment  (for  use  in  the  field,  or  in  a 
factory  environment)  which  can  measure  the  surface  stress  in  seconds.  An 
example  of  this  equipment  is  shown  in  Fig.  4^.  It  is  even  possible  to  use  only 
one  incident  beam  direction  and  examine,  simultaneously  with  two  detectors,  two 
points  on  the  diffraction  cone  which  come  from  grains  tilted  at  different  ^ 
values^.  However,  it  has  been  shown ^ that  in  this  case  the  stress  constant, 
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K,  is  much  larger,  which  implies  a larger  error  in  stress  from  a given  uncer- 
tainty in  29.  (Film  equipment  employing  this  "single  exposure"  technique  is 
still  sometimes  used  for  very  large  pieces.) 


Because  of  the  simple  form  of  Eqn.  6a,  it  is  particularly  easy  to  measure 
K by  applying  a series  of  known  elastic  loads  to  a specimen.  Tensile  or 
bending  loads  are  commonly  employed ^ . 


4.  The  General  Stress  State 

Because  the  measured  stresses  are  averages  over  a finite  depth  of  beam 
penetration,  we  should  consider  the  general  stress  tensor: 


*11 

*12 

*13 

*12 

*22 

*23 

*13 

*23 

*33 

Referring  to  Eqn.  1,  the  presence  of  e^-j  and  623  leads  to  curvature  in  the  d vs 
2 

sin  ijr  relationship,  as  their  contribution  depends  on  sin2^.  Such  terms  can  be 

important  in  machining,  or  in  wear;  that  is,  when  there  are  strong  components 

of  the  deforming  forces  tangential  to  their  surface,  and  strong  gradients  (since 

the  strains  e^>  ^23’  e33  must  vani-sh  at  the  surface).  Examples  of  this  effect 

(from  ref.  11)  are  shown  in  Figs.  5 and  6.  The  principal  effect  of  the  shear 

2 

stresses  is  to  cause  "d"  vs  sin  | to  be  non-linear  in  opposite  senses  for  positive 

and  negative  tilts.  This  ^ splitting  and  the  large  curvature  near  i(f“0  are 

distinguishing  features  of  this  phenomenon.  The  actual  stress  tensors  are 

given  in  the  captions  to  the  figures,  and  the  small  value  of  the  shear  stress 

to  cause  the  "1^  splitting"  is  particularly  noteworthy.  Also,  the  curves  for 

positive  and  negative  \|r  (or  by  changing  $ by  180°)  reverse  when  the  sign  of  the 

shear  stress  reverses , Fig.  6.  This  curvature  and  its  cause  were  first  reported 

(12) 

by  Walburger  , and  a simple  method  of  analysis  of  the  data  was  suggested  by 


and  a simple  method  of  analysis  of  the  data  was  suggested  by 
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Dolle  and  Hauk'-  . It  was  employed  in  an  extensive  study  for  the  first  time 
by  Dolle  and  Cohen , introducing  the  terms  a^  and  a^,  employing  Eqn.  1 
and  using  carats  to  indicate  averages  over  the  depth  of  penetration: 

al  = + e^J  ■ <e33>  + [<eu>cos2$  + <e12)sin2$  + <e22>sin2§  - <e33> ]sin2^ , 

(8a) 

a2  = " «#^.]  “ [<e13>cos$  + <e23)sin$]  sin|2*|.  (8b) 

2 

Thus,  (e33)  can  be  obtained  from  the  intercept  of  a^  vs  sin  this  value 

is  independent  of  §,  and  so  it  can  be  verified  by  measurements  at  different  § 

2 

tilts.  The  components  e^2>  ^22  can  136  obtained  from^a^/g  sin  At  $ = 

0°,  e33  " e33»  and  hence  e^,  is  obtained,  whereas  for  $ * 90°,  e22  can  136 

measured.  The  Sl.  component  is  then  determined  from  f.  ^ al 1 . From 

aa2  0 |sin2|  §-45 

in | 2,^1  e13  is  calculated  for  $ » 0 » «23  for  § » 90  . The  curves  in  Figs.  5 

and  6 are  the  fit  to  the  data  points  with  Eqn.  1. 

The  stresses  are  then  obtained  from  Eqn.  2. 

For  steels, 

measurements  generally  yield  a value  for  the  lattice 

parameter  of  2.8665(1)  X.  The  uncertainty  leads  to  a change  of  w + 30  MPa 
in  the  normal  stresses  and  hardly  any  variation  in  the  shear  stresses.  Measure- 
ments with  a well  annealed  powder  of  Cr  indicate  that  the  precision  is  15  MPa 
for  o32.>  ai2f  °22 * p33’  and  better  tban  8 MPa  for  o^,  Cf23»  Ibis  uncertainty 
is  often  less  than  the  variation  of  the  actual  stresses  from  point  to  point 
in  a specimen. 

Another  way  of  examining  the  stress  in  three  dimensions  is  with  high- 


In  this  procedure,  it  is^  necessary  to  know  "do." 


energy  X-rays  or  neutrons,  isolating  the  volume  element  to  be  examined  by 
suitably  narrowing  the  slits  defining  the  incident  and  diffracted  beams.  The 
feasibility  of  such  a procedure  is  considered  in  the  next  paper. 
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Except  for  materials  with  appreciable  content  of  a second  phase  (such 

as  a steel  with  more  than  « 0.4  wt.  pet  C) , some  kinds  of  deformation  may 

produce  strong  texture  at  modest  strains.  In  such  cases  oscillations  have 

2 

been  observed  in  d vs  sin  ^ as  shown  in  Fig.  7.  The  precise  interpretation  for 
this  effect  is  still  uncertain;  it  may  be  due  to  elastic  or  plastic  aniso- 
tropy^-*^. Tests  are  currently  underway  to  decide  which  factor  is  most  im- 
portant, after  which  procedures  may  be  possible  to  minimize  this  problem.  For  exam- 
ple, any  contributions  from  elastic  anisotropy  can  be  minimized  (for  any  materi- 


al) by  employing  hoo  or  hhh  reflections 


(5) 


5.  Gradients 

The  ability  to  detect  stresses  normal  to  the  surface,  as  was  discussed  in 
the  previous  section,  is  based  on  the  presence  of  gradients,  and  the  fact  that 
the  X-ray  beam  averages  over  its  penetration  (D) . Therefore,  it  is  appropriate 
to  examine^the  effects  of  such  gradients  because  .even  in  the  absence  of  stresses 
normal  to  the  surface,  there  will  be  effects  on  » (^22^*  Correc- 


’12' 


tions  are  available  for  this^’^',  and  also  for  stress  relief  due  to  layer 

removal if  the  gradients  are  examined  in  this  fashion.  But  the  presence  of 

2 

steep  gradients  can  also  be  detected  and  examined  from  the  form  of  "d"  vs  sin  $ 
without  layer  removal,  as  illustrated  in  Fig.  8.  For  this  figure,  the  stresses 


were  obtained  from  an  assumed  stress  profile  (in  the  z-direction) 


(5). 


<Oij)  “ aij  (z”°)  + j*  exp(-z/ T)gtj  (z)d  z 


’ij'  ij 

Here,  with  jj,  the  linear  absorption  coefficient: 

2 2 
sin  9 - sin  t 

T 2p,sin9  cos\jt  * 
for  ^ tilts  around  the  0 axis,  and: 

m sinQ  cos^ 


(9) 


(10a) 


(10b) 
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for  tilts  around  an  axis  parallel  to  the  plane  of  the  diffractometer.  The  term 

g„(z)  describes  the  variation  of  stress  with  depth  (z)  below  the  surface.  (Once 

these  average  stresses  are  evaluated  for  a particular  strain  distribution  (e^) 

2 

can  be  calculated  at  each  and  d vs.  sin  tjj  obtained  as  will  be  shown  below. 

Note  particularly  the  difference  in  curvature  near  ^ = 0 between  this  ef- 
fect (if  shear  stresses  (023)  * 0)  and  for  ^ splitting  (compare  Fig.  8 

with  Figs.  5 or  6).  Also,  the  values  of  d at  + * should  be  identical  in  this 
case,  so  the  effect  is  readily  distinguished  from  splitting. 

This  curvature  can  be  employed  to  obtain  information  about  gradients,  and 
we  now  present  a concrete  example  of  how  this  can  be  done,  with  data  from  a SAE 
1040  steel,  normalized  and  shot  peened.  We  make  the  following  assumptions  (all 
of  which  could  be  checked  in  a more  thorough  investigation) : 

a)  the  X-ray  elastic  constants  S,  and  S 2 do  not  vary  with  depth,  z: 

b)  varies  with  z only; 

c)  $ = 0°(by  definition) 

d)  cr,,(z)  " 022^^  (due  to  shot  peening  normal  to  the  surface); 


e)  f t j * 0 (no  appreciable  shear  deformation). 


Then: 


<633)  * ^S2(hk^)rsin2'Ko11)  + (1  - sin2^,)  (o33>]  + S1(hkX)[2<o11)  + <ff33>], 

$-0 


From  this  equation: 


<c33>*,*-0  " 2s1(hk^Xc11>  + <a33)  [i  S2(hkx)  + 8,0*0], 


. 2 

5 sin  ^ 


^S2(hkX)[(CTn>  - <533)], 


(12a) 


(12b) 


With  Eqns . 12,  first  estimate  of  and  (a^)  can  be  obtained,  fitting 

a straight  line  to  the  data.  Typical  data  is  shown  in  Fig.  9.  From  this  figure, 


I 


11 

these  first  estimates  are:  (c^)  “ “309  MPa,  (033)  ” 140  MPa.  Comparison 
of  Fig.  9 with  Fig.  8 assuming  the  gradients  are  linear  indicates  that  either 
g^  < 0 or  g^^  > 0,  or  „oth. 

The  measured  macroscopic  gradient  is  illustrated  in  Fig.  10a.  "t  can  be 

seen  that  the  stresses  vanish  at  «s  500  A simple  representation  of  this 

figure  was  taken  to  be  Fig.  10b.  The  terms  and  g„  were  varied  to  give  the 

2 

best  fit  to  the  data  : (c^),  {033)  > and  the  measured  slopes  of  "d"  vs  sin  ^ 

at  various  angles. 

From  the  average  (033)  * 140  MPa  and  the  penetration  depth,  j = 10 

2 

g^^  was  the  order  of  15  MPa.  From  the  curvature  of  "d"  vs  sin  it  could  be 

concluded  that  g^  was  -10  to  -15  MPa/^m.  These  values  were  actually  obtained 

simply  by  choosing  various  z ^ values  and  surface  stresses.  Of  course,  more 

sophisticated  least-squares  approaches  are  possible.  This  procedure  might  be 

particularly  fruitful  when  many  similar  samples  have  to  be  examined  for  dieir 

gradients,  and  the  gradient  in  one  specimen  has  been  measured  by  layer  removal. 

Finally,  it  is  worth  noting  that  peak  shapes,  as  well  as  positions  have 

(14) 

been  employed  by  Murakami v ' to  obtain  information  on  gradients  in  thin  foils. 
While  his  assumptions  as  to  the  stress  tensor  are  too  restrictive  for  general 
use,  his  procedure  could  be  expanded. 

6.  Factors  Affecting  Accuracy 

We  have  already  mentioned  one  possible  source  of  inaccuracy,  the  elastic 
constants.  The  actual  error  in  stress  from  the  effects  to  be  discussed  below 
can  be  evaluated  through  the  stress  constant  K and  Eqns.  6,  or  by  employing 
experimental  X-ray  elastic  constants  and  S^.  A particular  error  in  29, 
when  multiplied  by  K for  the  chosen  material,  gives  a first  estimate  of  the 


error  in  stress. 


12 


L 


All  data  should  be  corrected  for  the  variation  in  scattering  factor  and 

o 


Debye-Waller  factor  if  the  peak  is  very  broad  (>  « 5 in  half  breadth). 
Corrections  for  the  Lorentz  polarization  and  absorption  factors  (LPA)  can  be 
written  for  each  position,  29,  across  a peak.  Assuming  filtered  radiation: 


LPA 


1 + cos  29 


• 2„ 

sin  9 


(1  - tan\|j  cot9). 


(13a) 


for  tilts  around  the  29  axis,  and: 


LPA 


1 + cos  29 


sin29 


(13b) 


for  tilts  around  an  axis  parallel  to  the  goniometer.  The  latter  method 
is  increasingly  employed  in  Europe.  It  is  not  limited  to  |ty|  £ 9,  as  is  a 
tilt  around  the  9 axis.  Furthermore,  the  path  lengths  of  the  incident  and 
diffracted  beams  within  the  specimen  are  equal,  so  absorption  corrections  are 
not  involved  (Eqn.  13b  vs  Eqn  13a).  It  is  a helpful  technique  for  measurements 
in  difficult  locations  such  as  on  gear  teeth,  where  the  incident  or  diffracted 
beam  can  be  blocked  by  the  gear  itself. 

Because  these  corrections  are  small  and  straightforward  there  is  no  need  to 
avoid  them,  and  therefore  no  need  to  estimate  errors  from  these  sources. 

Locating  the  specimen  over  the  center  of  rotation  can  be  particularly 
troublesome  especially  when  the  geometry  of  the  specimen  is  complicated  and  this  is 
one  of  the  largest  sources  of  inaccuracy.  There  are  two  ways  of  achieving  this 


location.  The  lattice  parameter  for  two  or  more  different  reflections  (widely 

:d  in  2( 

(17) 

parameter'  : 


cos^Q  o 

separated  in  29)  can  be  plotted  vs  — r-.  Then  (for  i|(»0  ),  with  "a"  the  lattice 

s lriu 


-21  cos  9 > 
a R sin9 


(14) 


with  5 the  sample  displacement  and  R the  goniometer  radius.  The  displacement 
can  therefore  be  calculated  from  the  slope  of  such  a plot,  and  the  sample 
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displaced  appropriately. 

Another  way,  particularly  important  when  ^ splitting  is  suspected, 
is  to  examine  the  peak  position  of  a stress  free  powder  (painted  in  a thin 
layer  on  the  specimen  in  an  acetone  solution)  at  a fixed  29  (close  to  the 
reflection  angle  of  the  specimen),  and  then  varying  y . An  example  of  such  a method 
is  shown  in  Fig.  11. 

/ iq\ 

In  Japan,  a parallel  beam  method  is  employed''  J , to  eliminate  the  sensi- 
tivity to  sample  position.  In  Fig.  12,  a comparison  of  the  peak  shift  due  to 
sample  displacement  is  made  between  this  method,  the  parafocussing  method  (for 
which,  the  receiving  slits  are  moved  to  the  correct  focus  at  each  , and  a 
stationary  slit,  divergent  beam.  The  parallel  beam  technique  is  obviously 
least  sensitive.  However,  the  peak  is  broadest  in  this  case  and  the  time 
for  data  analysis  is  sa  1/3  longer  for  the  same  statistical  error  (Section  7) 

(19) 

than  parafocussing,  and  the  reproducibility  twice  as  poor  for  sharp  profiles 
The  stationary  slit  method  appears  to  be  the  best  compromise.  However  in  this 
case  (but  not  in  parafocussing)  the  shape  of  the  peak  may  vary  with  ^ tilt. 

How  this  alters  various  methods  for  peak  location  is  illustrated  in  Fig.  13. 

Bias  corrections  may  be  taken  from  this  plot. 

(If  parafocussing  is  employed,  and  the  sample  has  a strong  texture,  it 
is  important  to  introduce  a permanently  fixed  slit, defining  the  beam, along  the 
9 axis;  it  is  not  moved  when  the  scanning  slit  is  moved  to  the  correct  focus. 

This  assures  that  the  same  portion  of  the  diffraction  cone  is  seen  at  all 
tilts(20).) 

An  excellent  presentation  of  the  bias  from  various  instrumental  misalign- 
ments can  be  found  in  ref.  2;  while  the  equations  in  this  reference  are  based 
on  the  center  of  gravity  of  a peak,  they  can  be  used  as  a guide  for  any  definition 
of  the  peak.  Another  source  of  analytical  approximations  to  the  errors  is  ref.  22. 

L ^ 
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For  this  paper,  we  have  made  computer  calculations  of  the  peak  shifts  by  a ray-tracing 

method.  The  intensity  of  the  beam  was  assumed  to  be  uniform  across  its 

cross  section.  This  over-emphasizes  the  error  and  the  resultant  shifts  were 

multiplied  by  1/3,  in  agreement  with  the  findings  in  similar  comparisons 

(13) 

by  Zantopulos  and  Jatczak  . The  errors  summarized  in  Table  I are  the 
shift  in  29  for  a ^ tilt  from  +45°  to  -45°.  If  only  +i|f  is  employed,  the 
errors  are  reduced  by  an  additional  factor  of  1/2.  This  table  can  serve  as 
a quick  guide  to  the  important  aspects  in  alignment  of  the  diffractometer 
for  stress  determinations.  Table  II  gives  some  typical  errors  for  the  two- 
tilt  method  and  parafocussing,  with  tilt  around  the  0 axis.  Clearly,  the 
bias  due  to  geometric  errors  can  be  held  to  quite  small  values,  with  good 
alignment  and  a judicious  choice  of  slits. 

7.  Factors  Affecting  Precision 

Automation  of  stress  measurements  has  been  achieved  in  several  locations. 

One  such  program,  for  example,  includes  sample  alignment  over  the  tilt  axis,  and 

preliminary  scans  to  locate  the  peak  and  to  determine  the  time  to  achieve 

(19  24) 

an  operator-specified  precision  ’ . Such  on-line  controls  allow  studies 

to  be  made  of  factors  affecting  reproducibility. 

As  shown  in  Table  II  the  error  in  determining  a peak's  location  is  one  of 
the  largest  sources  of  uncertainty  in  the  results.  A number  of  methods  have 
been  presented  to  define  the  angle  of  diffraction,  29,  from  the  diffraction 
profile,  including  the  apex  of  a parabola  or  cubic  polynomial  fit  to  the  upper 
portion  of  the  profile^  , the  centroid v J and  the  half  value  breadth'’  . 


Comparisons  between  these  methods  indicate  that  the  apex  cf  a least  squares 
parabola  fit  to  that  portion  of  the  profile  above  8571 i of  cne  maximum  intensity 
usually  provides  the  most  reproducible  procedure  for  defining  the  diffraction 
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(19  24) 

profile  when  determining  peak  shifts  ’ . This  is  shown  in  Table  III. 

If  the  intensity  1^  accumulated  at  each  jth  interval  in  29  is  taken 
at  2n+l  observation  points  (the  center  point  being  defined  as  the  working 
origin  29q)  in  equal  increments,  5,  of  29,  the  peak  location  in  a parabolic 
fit  is  given  by^^l 


where : 


28  - 28o  - (»22-non4)62M1/2r.2(n262M<j-»<)M2), 


n „ i+1  n . i 

n.  - E j andM.  M Z J 1 = 

j“-n  j**-n  J 


(15a) 


(15b) 


(Angular  intensity  corrections  such  as  the  Lorentz  polarization  and  absorp- 
tion factors  must  be  performed  on  the  raw  data  before  determining  the  peak 
position.) 

Scatter  in  the  peak  position  arises  from  statistical  fluctuations  in 
the  X-ray  intensities  and  the  precision  of  the  measured  value  is  properly 
evaluated  by  its  standard  deviation  due  to  counting  statistics.  The  standard 
deviation  in  peak  location  (Eqn.  15)  due  to  the  random  counting  errors 
is  (19, 24); 

4 2 2 

r 6 (n2~n  n, ) n 3 2 2 2 

S(2V  “ ITT  — -4  2 (V  *1*  + (V  VV96  J 

4n2(n26  MQ-nM2)  jm_n 


-n^3^)2  (16) 


where  the  variance  in  the  power,  S (1^)  depends  on  the  method  of  data  accu- 
mulation and  is  given  by: 


SPT<V  • li/C 

4c<v  - 


for  fixed  time,  t, 


for  fixed  counts,  c. 


(17a) 


(17b) 


The  error,  S(29  ),  introduces  an  error  in  the  calculated  interplanar  spacing 
P 

given  by: 
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[I 


S(d) 


j-XcosQ  - S^29p^ 

2,inV  * 


(18(P* 


(18) 


where  X represents  the  characteristic  wavelength  and  S(29)  is  in  degrees  29. 
The  subsequent  statistical  counting  error  in  measured  residual  stress 


is: 


Z (sin2’jr-sin2'j[)2» S2(di|r) 

t 

r . . 2 -r-2  .2.2 

[E  (sin  ii-sin  ■(,)  ] 


[vJS"].  <i9> 


for  the  sin  ^ technique  (Eqn.  5)  and: 


S(o$)  - K[S2(29q)  + S2(29+)3^, 


(20) 


for  the  two  tilt  method  (Eqn.  6). 

Fig.  14a  shows  the  distribution  of  errors  in  a large  number  of  actual 

measurements  when  0.01°29  was  chosen  as  the  desired  statistical  precision 

for  each  peak  in  a control  program  employing  the  above  equations.  The  half- 

-3  o 

width  of  the  distribution  is  only  2.4  *10  29.  But  when  the  distribution 

of  errors  in  repeated  measurements  on  the  same  specimens  was  examined  (without 

moving  it),  the  results  in  Fig.  14b  were  obtained.  Clearly,  the  scatter 

(19  24) 

is  greater  than  the  statistical  error.  In  a series  of  experiments  * , 

the  factors  required  to  achieve  reproducibility  equivalent  to  the  statistical 
error  were  examined.  The  results  can  be  summarized  as  follows: 

a)  Unless  the  peak  is  very  sharp,  at  least  7 points  should  be  employed 
in  a parabolic  fit,  if  the  two-tilt  method  is  employed  - not  the  three  points 
usually  taken.  (The  total  measuring  time  can  be  kept  fixed.) 

b)  If  the  peak-to-background  ratio  is  low  (2:1)  or  the  peak  is  broad 
(2°  or  more),  background  subtraction  is  required,  at  least  in  the  first  scan 


to  determine  the  region  of  fit.  (This  was  not  done  for  the  data  in  Fig.  14 
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and  this  is  the  main  reason  for  the  greater  than  expected  precision.  ) 

Even  with  (a)  and  (b)  the  scatter  may  be  twice  that  predicted  by  the 

above  statistical  analysis  for  the  two-tilt  procedure.  But: 

2 

c)  With  the  total  time  kept  fixed,  the  sin  ^ method,  with  at  least  a 
three  point  parabolic  fit  and  background  subtraction  will  achieve  the  expected 
error. 

d)  An  alternative  to  (c)  is  to  employ  a position  sensitive  detector 
(PSD)  and  30  or  more  points  in  the  region  of  fit.  In  this  case,  background 
subtraction  is  not  necessary. 

Automation  alone  will  reduce  the  time  of  measurement  by  2.  A PSD 

(24  27) 

will  result  in  time  saving  of  a factor  of  5 or  more  * . Indeed  if  an  error  of 

34  MPa  (5  ksi)  is  acceptable,  it  has  been  shown  that  in  some  cases  measurements 
can  be  conq>leted  with  the  two  tilt  method  with  2 seconds  for  each  tilt^\ 

Finally  it  is  instructive  to  consider  the  total  error  possible  in  a 
residual  stress  measurement  due  to  all  the  factors  discussed.  Samples  cir- 
culated to  many  investigators  can  give  perhaps  the  best  view  of  the  combined 
accuracy  and  precision  to  be  expected.  The  following  data  is  taken  from  ref. 

28.  The  two  tilt  method  was  employed,  with  steel  specimens. 

1.  SAE  Round  Robin  No.  3 (11  laboratories  reported,  1959) 

a.  flat  specimen  (broad  peak) 
stress  level:  14  MPa 

standard  deviation  among  laboratories:  + 10.3  MPa 

b.  1010  annealed  flat  steel  (sharp  peak) 

.tress  level:  -5  MPa 


standard  deviation:  + 16.5  MPa 


18 


Ik 


2.  SAE  Round  Robin  No.  4 (25  laboratories,  including  many  with  little  or 
no  experience  in  residual  stress  measurements) 

a.  flat,  shot-peened  (R£  63) 
stress  level:  593  MPa 
standard  deviation:  + 41  MPa 

b.  1045  Round  bar,  6.22  cm  diameter  (R£  61/62) 
stress  level:  910  MPa 

standard  deviation  (axial  direction):  + 56.5  MPa 
standard  deviation  (longitudinal  direction):  + 72  MPa 
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Table  I 

The  Influence  of  Various  Mis-Alignments  on  the 

Accuracy  of  Stress  Analysis 

o 

+ - large  error  (0.07 

or  more) 

o 

- ■ small  error  (0.01 

1 

or  less) 

m m medium  error 

f tilt 

* tilt 
around 

around 

axis  in 

type  of  error 

9 axis 

diffraction 

plane 

centers  of  0,  ijr  rotation 
not  coincident  (a  » 0.1  mm) 

o 

4 horizontal  divergence 
o 

4 vertical  divergence 

1°  horizontal  tilt  of 
beam  to  the  zero 
direction 


1 vertical  tilt  of  beam 
to  true  zero  direction 

tube  mis-alignment 
(Ax,  Ay,  Az  1 mm) 

displacement  of  sample 
from  the  tilt  axis 
(0.2  mm) 

sample  curvature  (beam 
much  smaller) 


m to  + 
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Table  II 

Typical  Bias  in  Peak  Location  for  20  * 156° 
(Approximately  the  211  CrK  Peak  from  Iron) 


peak  shift  between  peak  shift  between 

cause  >(r  “ 0°  and  * 45  =*  0 and  ” 60 


°20 

MPa 

°20 

MPa 

peak  location  + O.O1°20 

+ .02 

11.9 

+ .02 

+1700 

total  horizontal  beam 
divergence  of  1 

-.0006° 

-0.35 

-.0025 

-215 

vertical  beam  divergence ^ ^ 
(assuming  strong  texture 
using  divergent  Soller  slit; 
no  receiving  Soller  slit) 

+.002 

+1.20 

+ .002 

+170 

sample  displacement, 

AX’  * +.025  ram 

+. 0034 

+2.1 

+.0088 

+760 

i^-axis  displacement, 

AX'  ■ +.025  mm 

+ .002 

+1.2 

+.0068 

+585 

, *+• 

maximum  total  errors 

a)  in  -20  direction 

b)  in  +20  direction 

-.008 

+.0068 

-4.8 

+4.1 

-.0201 

+.0149 

-1150 

+860 

Note:  Maximum  error  is  either  one  of  these  but  not  the  total  range. 

+ Calculated  for  steel  from  a = K (20-20  ) where  = 602  MPa/°20  and 
K^q“40B.1  MPa  ksi/°20.  Does*not  ^include^error  in  peak  location  as  it  is 
dependent  on  time  of  data  collection. 


Table  III 

Precision  of  Various  Measures  of  Profile  Position 
(10  measurements) 


Sample 

Time 

(sec) 

FWHM* 

(°26) 

Half-Value  Breadth 
(°20) 

Centroid 

(°26) 

Parabola 

(°20) 

AISI 

1090-1 

50 

.45 

156.149  (+  .021) 

156.096  (+  .011) 

156.186  (+  .016) 

AISI 

1045-2 

100 

3.45 

155.336  (+  .064) 

155.396  (+  .085) 

155.413  (+  .020) 

FWHM  is  the  full  width  at  half  of  the  maximum  intensity 

The  term  in  ( ) represents  one  standard  deviation  from  the  average  position  over  the 
10  measurements. 
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Fig.  1 


Fig.  2 

Fig.  3 


Fig.  4 


Fig.  5 


Figure  Captions 

Top:  Schematic  of  a diffractometer.  The  incident  beam  diffracts 

from  planes  that  satisfy  Bragg's  law,  in  grains  with  these  planes  paral- 
lel to  the  sample's  surface.  If  regions  near  the  surface  are  in  com- 
pression, with  the  stress  parallel  to  the  surface,  because  of  Poisson's 
ratio,  these  planes  are  further  apart  than  in  the  stress-free  state. 

The  "d"  spacing  is  obtained  from  the  peak  in  intensity  vs.  scattering 
angle  20  (insert),  and  Bragg's  law. 

Bottom:  After  the  specimen  is  tilted,  diffraction  occurs  from  other 
grains,  but  from  the  same  planes,  and  these  are  less  separated  than  in 
(a),  because  there  is  a component  of  the  residual  stress  normal  to  the 
planes. 

The  axial  system.  The  are  the  laboratory  system  and  measurements 

of  Ad/d  are  made  along  L^.  The  describe  the  specimen. 

2 

Interplanar  spacing,  "d",  vs  sin  ty,  normalized  AISI-1010  steel, 

deformed  in  tension  along  to  a true  strain  of  8 pet.  211  peak, 

CrK  radiation;  « .o * -107.1  MPa  (-15,303  psi).  From  ref.  22. 
a $=u 

The  PARS  system^  - a portable  (X-ray)  analyzer  for  residual  stresses. 
A)  linear  position  sensitive  detector  to  record  the  entire  diffraction 
profile  without  detector  motion;  B)  air-cooled  miniature  X-ray  tube 
(50  kv,  2 mA) ; C)  exit  soller  slits;  D)  rods  to  position  unit  with 
respect  to  a specimen;  E)  one  (of  two)  handles. 

Normalized  plain  carbon  steel  strip  (0.6  wt  pet  carbon)  ground  under 
flowing  water  to  remove  5 per  pass^^.  Various  $ tilts  shown. 


The  curves  are  fits  to  the  data  points  with  eqn.  (1).  The  stress 
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Fig.  6 


Fig.  7 


tensor  in  MPa  is 

^ 390  14 

14  306 

63  -1 

211  peak,  CrK  radiation. 
a 

2 

(e^)  vs  sin  f.  Normalized  plain  carbon  steel  strip  (0.6  wt  pet  carbon) 

ground  to  remove  5 p,m  per  pass  under  flowing  water § = 0°.  Strain 

-3 

in  units  of  10  . Specimen  C2  was  ground  only  in  the  direction  of  the 
+P^  axis.  For  specimen  C5,  the  direction  of  grinding  was  reversed  each 
pass  with  the  final  pass  opposite  to  that  for  C2  (in  the  -P^  direction). 
The  stress  tensor  for  C5  is  given  in  Fig.  5.  For  C2  (in  MPa)  it  is: 


*3  \ 


/ 199  -10  -63  \ 

-10  86  5 

^ -63  5 84 

Note  that  the  curvature  for  y > 0,  | < 0 reverses  when  reverses. 

The  curves  are  the  fits  to  the  data  points  with  the  stress  tensors  and 

Eqn.  1.  0:  < 0;  • i)r  > 0. 

211  peak,  CrK  radiation 

a 

2 

Interplanar  spacing,  "d",  vs  sin  | for  a CU3AU,  specimen  pulled  in 
tension  along  P^^  to  a true  strain  of  35  pct^22\  420  peak,  CuK  radiation. 
The  term  f(0,t)  is  the  area  under  the  Bragg  peak  at  each  \|i  position. 


2 

Fig.  8 The  effect  of  various  gradients  on  "d"  vs  sin  ty.  The  gradient  is 

assumed  to  be  linear  over  the  X-ray  beam's  depth  of  penetration,  with 
slope  . Typical  values  were  assumed  (for  0^,  200-600  MPa  over  20  ^.m, 
for  033,  100  MPa  per  um,  (zero  at  the  surface)). 
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Fig.  9 


Fig.  10 


Fig.  11 


Fig.  12 


Fig.  13 


Fig.  14 


2 

Measured  'fl"  vs  sin  y,  shot  peened  (normalized)  AISI  1040  steel;  data 
averaged  over  nine  specimens.  Error  bars  from  observed  variation  in 
29  in  repeated  measurements.  CoK^,  310  reflection. 

a)  The  gradient  measured  by  layer  removal  for  one  of  the  specimens  in 
Fig.  9 (Etchant:  100  parts  30  pet  1^0,  10  parts  48  pet  HF). 

Data  uncorrected  for  layer  removal  and  depth  of  penetration.  310  peak,  CoK  radi 

b)  Assumed  gradients  for  calculations  described  in  text  to  fit  data  in  (a) . 

2 

The  effect  of  sample  displacement  on  "d"  vs  sin  . Measurements  made 

with  an  annealed  Cr  powder  painted  on  a steel  specimen  in  an  acetone  solution. 

211  reflection,  CrK  radiation. 

a 

Measured  effects  of  sample  displacement,  for  various  focussing  techniques 

employed  in  stress  measurements.  Normalized  1045  steel,  211  reflection, 

CrK  radiation. 

a 

Calculated  peak  shift  due  to  K - doublet  for  various  definitions  of  the  pea! 

*1  2 

a)  Special  slits  to  symmetrize  peak  (Wolfstieg,  U.,  Die  Symmetrisierung 
unsymmet rischer  Interferenzlinien  mit  Hilfevon,  Spezial  Henden, 

Harterei-Techn.  Mit.  31^,  23-26,  (1978)). 

b)  center  of  gravity 

c)  middle  of  half-width 

d)  parabolic  fit  to  top  15  pet 

The  same  Gaussian  functions  were  assumed  for  the  a 2 components. 

211  reflection,  CrK  radiation. 

o 

a)  Precision  of  peak  location  (parabolic  fit,  top  15  pet).  572  peaks  ; 
actual  standard  deviations  when  the  indicated  precision  was  requested. 

b)  440  peaks  each  examined  4 times;  indicating  the  spread  in  each  of 
these  4 repeated  measurements. 

Armco  iron  and  steels,  211  reflection,  CrK  radiation. 

a 

7 point  parabolic  fit  for  peak  positions. 
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